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Abstract. We present the EROS nearby supernova (z ~ 
0.02 — 0.2) search and the analysis of the first year of 
data (1997). A total of 80 square degrees were surveyed. 
Eight supernovas were detected, four of which were spec- 
troscopically identified as type la supernova?. The search 
efficiency was determined with a Monte-Carlo simulation 
taking into account the efficiencies for both supernova de- 
tection and host galaxy identification. Assuming that for 
a given galaxy the supernova rate is proportional to the 
galactic luminosity, we compute a type la supernova ex- 
plosion rate of: K = 0.44t°; 35 h 2 /10 10 L 0B / 100 yrs 
at an average redshift of ~ 0.1 where the errors are re- 
spectively statistical and systematic (type misidentifica- 
tion included). 

Key words: supernovae, rates, cosmology 



1. Introduction 

Type la supernovae have been shown to be powerful cos- 
mological distance indicators useful for the determination 
of the expansion rate Hq and the density paramete rs for 
matter and vacuu m ene rgy fijf and Qa (Riess et al. 1998 : 
Perlmutter et al. 1999 ). As such, they are the subject of 
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intense study and a number of nearby and distant su- 
pernova searches have been launched. Studies of type la 
light-curve shapes, peak magnitudes, colors and spectra 
have been performed in order to get better insight into 
the nature of these exploding stars. 

Supernova rates, and in particular their evolution as 
a function of redshift, is another probe provided by su- 
pernova science. Not only are these rates important for 
tracking the chemical evolution of the universe, but as a 
tracer of stellar evolution, they also contain information 
on which type of progenitor system produces type la su- 
pernovae. If the evolution of the supernova rate is under- 
stood, they may also be used to determine the cosmologi- 
cal parameters Qm arL d Q\ via the number count-redshift 
relationship. With enough statistics, it is possible to mea- 
sure the supernova luminosity function which is of critical 
importance for understanding Malmquist bias. 

The supernova rate has been measured at low red- 
shift in searches based on visual scanning of photographic 
plates (Cappellaro et al. 1997). At high redshift the rate 



has been measured with automatic subtraction of CCD 
images (Pain et al. 1996). The supernova rate for nearby 



galaxies has previously been estimated using CCD data by 



Muller et al. (1992). Here we report the first measurement 



of the rate for low redshift galaxies using CCDs and with 
a full numerical calculation of the detection efficiency. The 
search was performed by EROS (Experience de Recherche 
d'Objcts Sombres). 
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In 1990 EROS engaged in a search for massive com- 
pact halo objects (MACHOs) via gravitational microlens- 
ing. Since 1996 EROS has been using a dedicated 1 me- 
ter telescope, the Marly telescope, at the ESO La Silla 
observatory. The telescope is equipped with two detector 
mosaics, each consisting of eight 2048x2048 pixel CCDs 
covering a 0.7 x 1.4 deg 2 field. A dichroic cube allows the 
simultaneous recording of two images in wide (AA ~ 200 
nm) non-standard red (A ~ 760 nm) and visible (A ~ 560 
nm) bands hereafter referred to as i?Eros and Ve T os- The 
telescope, came ra an d telescope operation are described 
in Bauer et al. (|l997j ). 

The greater part of EROS observing time is devoted to 
the search for microlensing events by observing the Mag- 
ellanic Clouds, the Galactic Bulge an d the Galactic Disk 
(see e.g. Palanque-Delabrouille et al. 1998| ; Derue et al. 
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1999j) . However, EROS also devotes 10% of its observing 
time to the observation of fields at high Galactic latitudes 
in order to search for supernovae at low redshifts and to 
measure stellar proper motions. 

The EROS supernova search aims at discovering a ho- 
mogeneous sample of supernovae in the low redshift range. 
Our main scientific goals are the derivation of the nearby 
supernova rates, and the study of the empirical correlation 
of the peak luminosity with the light-curve shape. EROS 
has discovered more than 60 supernovae since 1997. 

In this paper, we present a measurement of the type 
la supernova explosion rate. Because of backgrounds due 
to variable stars and asteroids, we restricted the search to 
identified galaxies. To derive a total supernova rate, we 
therefore assume that the rate for a given galaxy is pro- 
portional to the galactic luminosity. For this low statistic 
supernova sample, no attempt is made to derive rates for 
different galaxy types. 

In Section 2, we give the details of our supernova search 
that are relevant for the rate determination. In Section 3 
we present the subset of data considered in this analysis 
and the supernovae we found. The galaxy sample for the 
supernova search is described in Section 4. In Section 5 
we describe the determination of our detection efficiency 
using a Monte Carlo simulation and derive the supernova 
rate. Finally, in section 6 we present a discussion and inter- 
pretation of our results. More details on this measurement 
can be found in Hardin ( |1998| ). 



2. The search technique 

We searched for supernovae by comparing an image of a 
given field with a reference image of the same field taken 
at least two weeks before. The images are taken during 
dark time. Only the visible band images are used for the 
search, but we use both bands when computing the abso- 
lute calibration of the images. The technique is illustrated 
in Fig. [I]. Each observation consists of two five minute ex- 
posures which are combined to form a ten minute exposure 
after identification of cosmic rays. Five to ten fields were 



Fig. 1: Illustration of the detection technique (detection 
images of SN 1997eb). a) Reference image taken on Oc- 
tober 7, 1997. b) Detection image taken on November 7, 
1997. The supernova can barely be distinguished from the 
galaxy nucleus, c) The subtracted frame, obtained after 
geometric and photometric alignment of a) and b) frames. 
The stellar PSFs are matched by convolving the best im- 
age. The supernova now stands alone. 



observed each night. For the observations used in this pa- 
per, calibration fields (Landolt 1992) were also observed 
before and after the series of supernova fields. This allowed 
us to calibrate our non-standard filters and to eliminate 
non- photometric nights from the rate measurement. 

The automated processing permits the remote analysis 
in France on computers set up at the La Silla observatory 
of up to 20 deg 2 during the day following the night of ob- 
servation. Since the processing is monitored from France, 
it requires robustness and a limited volume of output. 

The method for comparing the reference and the search 
images is as follows. The reference image is aligned geo- 
metrically and photometrically with the search image, and 
the image of superior seeing is convolved to match the 
point spread function (PSF) of the other image. The two 
images are then subtracted and the frame thus obtained is 
searched for star-like sources. The selection of the candi- 
date supernovae is performed by applying a set of cuts on 
the objects detected on the subtracted frame. These cuts 
are tuned with a Monte-Carlo simulation. 

The first and most important cut is applied on the to- 
tal flux of the candidate, and is equivalent to a limiting 
magnitude of Vj ~ 21.5. Afterwards, a series of cuts in- 
volving the comparison of the candidate shape (i.e. the 
second moments of a Gaussian fit) with that of the frame 
stars is designed to eliminate inaccurate subtractions. Fi- 
nally, the detection threshold is set slightly higher within 
1.8 arcsec of the center of bright galaxies. 

The treatment ends with the visual scanning of ~ 5 
candidates per square degree by two independent ob- 
servers. Most of the candidates are residual cosmic rays, 
subtraction artifacts, asteroids, and variable stars. The 
residual cosmic rays are identified by the fact that they 
only appear on one of the two five minute exposures. As- 
teroids and variable stars are eliminated by requiring the 
candidate to appear in or near a host identified as a galaxy. 
As a consequence of this last requirement, no supernovae 
with an undetected host can be found. This is taken into 
account in the calculation of the search efficiency. 
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3. The data set 

For this analysis, we have used the data obtained dur- 
ing the two search runs conducted during October and 
November 1997. A description of the observed fields is 
given in Tab. [I]. They have been divided into 3 zones. The 
southern hemisphere fields have been chosen far from the 
Galactic plane (mean Galactic latitude <6>^ —70°), in 
a region covered by The Las C ampanas Redshift Survey 
(LCRS, Shectman et al. |l996|) . While the Abell Cluster 
fields^ The northern hemisphere fields have been chosen 
so as to be observable from the La Silla (Chile) and Apache 
Point (New Mexico, USA) observatories at the time of the 
November 1997 run. As a consequence, their mean Galac- 
tic latitude (< b >^ —50°) is not as high as in the southern 
hemisphere zone. The total area of sky (80 square degrees) 
was covered over a period of 10 days. 

Among the 8 detected supernovse, 5 could be classified 
according to their spectra. The type of the three remaining 
supernovae was obtained using the shapes of their light 
curves. The characteristics of these eight supernovae are 
summarized in Tab. 2. 



4. The galaxy sample 

In order to determine the supernova rate, we must have 
a well-defined sample of potential host galaxies with a 
known distribution of redshifts. In this section, we describe 
our galaxy identification and redshift determination algo- 
rithms and their verification using LCRS galaxies in our 
fields. 



4-1- Galaxy identification 

The position and flux of the sources on the reference 
and detection frames are measured with the software for 
source extraction SExtractor (Bcrtin & Arnouts 1996| ). 
SExtractor measures the total flux of extended sources, 
and provides in addition a star-galaxy separation estima- 
tor, CLASS STAR, which uses the intensity profile. 

Our galaxy selection criterion is based on the (total) 
Vferos magnitude and CLASS star. Because of the PSF 
variation over the CCD mosaic, the cut applied on CLASS 
STAR differs from one CCD to the other. On the other 
hand, the cut on Ve TO s is set equal for all frames and 
corresponds to a magnitude R c0 ~ 18.7 for the October 
1997 search fields, and R c0 ~ 18.4 for the November 1997 
search fields. This difference is due to Galactic absorp- 
tion. The cut on the galaxy magnitude was to ensure that 
the host was classified as such during the visual scanning. 
At z = 0.1, this corresponds to an absolute magnitude 



EROS galaxy counts - comparison with published counts 




Our original program included a st udy o f the peculiar ve- 



locities of these clusters (Reiss et al., 1998). are centered on 
Abell Clusters at a mean redshift of z ~ 0.17, the search in 
these fields was not restricted to the galaxies that belong to 
these clusters. 



21 22 23 
R c magnitude 

Fig. 2: Comparison of EROS galaxy counts (histograms) 
with counts from Bertin & Dennefeld (1997; stars), Weir et 
al. (1995; bold dash-dot line), and the northern (bold full 
line) and southern counts (bold dashed line) from Picard 
( 199ip . The dotted line histogram shows all the detected 
sources on the frames, both stars and galaxies, the full 
line histogram the objects selected as galaxies (with R c < 
20.2), and the hatched histogram the objects selected as 
potential hosts for the simulation, i.e. with a magnitude 
R c ~ 18.7. This corresponds to a solid angle of 38 deg 2 
(October 1997 data). There is good agreement between 
our galaxy counts and published counts. 



of M Rc 20. + 51og(i/ /60kms- 1 Mpc _1 ) for the host 

galaxy. 

This selection criterion was first tested on the fields 
shared with the LCRS: 90% of the LCRS redshift catalog 
galaxies were identified as such by our criterion, and 95% 
of the selected objects with a magnitude R c < 17.5 (LCRS 
completeness limit) were LCRS galaxies. 

Having identified the galaxies on our frames and 
derived their R c (total) magnitude from their VE ros and 
i?Eros magnitudes and the calibration (see section 2), we 
obtain a galaxy count per magnitude and per deg 2 , which 
we can compare with published counts. The agreement is 
illustrated in Fig. ||. One can also identify the comple- 
mentary set of objects as stars and obtain in this way star 
counts as a function of magnitude. These star counts com- 
pare well with the Galactic star-count model of Bahcall & 
Soneira (|l984| ). 

At this stage, ~ 225 objects per deg 2 are selected as 
potential host galaxies. 
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zone 


# fields 


a (J2000) 


5 (J2000) 


I 


6 


October 1997 run 


Southern hem. 


43 


22hW <a< 2Zi30' 


8 ~ -40° 


250° < Z < 360° 


-75° < 6 < -60° 


Abell Cluster 


8 


2Zi00' <a< 3ft00' 


5-0° 


150° i I £ 180° 


-50° < 6 < -60° 


November 1997 run 


Abell Cluster 


4 


2h00' <a< 3ft00' 


5 ~ 0° 


150° < Z < 180° 


-50° < b < -60° 


Northern hem. 


36 


3Zi00' <a< 4ft00' 


5 ~ 0° 


175° S Z < 190° 


-35° < 6 < -55° 



Tabic 1: List of the surveyed sky zones during the October and the November 1997 runs. The number of fields surveyed 
in each zone is given in column 2. Due to a temporary malfunction of CCD jj= 7 on the camera corresponding to the 
visible band, one field only covers 0.871 deg 2 . The total surveyed area is thus ~ 45 (resp. ~ 35) square degrees for 
the October (resp. November) 1997 run. The southern hemisphere fields are included in the Las Campanas Redshift 
Survey fields. 



IAU name 


date 


type 


2 


host 
R c mag. 


discovery 
Vj mag 


dist. from 
host core 


SN 1997dh (1) 


20/10/1997 


Ic 


0.05 


15.7 


20. 


5.0 " 


SN 1997dk (2) 


26/10/1997 


la 


0.05 


15.2 


20. 


11.0 " 


SN 1997dl (2) 


26/10/1997 


la 


0.05 


15.2 


19.5 


11.0 " 


SN 1997dm (2) 


26/10/1997 


Hp (l.c.) 


0.03 


14.4 


21. 


6.7 " 


SN 1997eb (3) 


19/11/1997 


II (l.c.) 


0.08 


17.1 


21.5 


2.9 " 


SN 1997ec (3) 


20/11/1997 


iip (i.e.) 


0.12 


16.3 


21.5 


3.7 " 


SN 1997ed (3) 


22/11/1997 


Ia 


0.15 


17.5 


21. 


2.8 " 


SN 1997ee (3) 


23/11/1997 


la 


0.17 


17.0 


21. 


8.6 " 



Table 2: Main characteristics of the supernova? detected during the October and the November 1997 runs. In the third 
column, (l.c.) indicates that the supernova type was derived from its light curve. All magnitudes are corrected for 
Galactic absorption. Note that SN 1997dk and SN 1997dl occurred in the same galaxy 2 months apart according to 
their spectra. The references are as follows: (1) IAU6760, $997} (2) IAU6762, |i~997|; (3) IAU6782, |l997|, and IAU6785, 



1997 



4-2. Magnitude-redshift relationship 

For a galaxy of a given magnitude -R C gai, the probability 
that its redshift is z is given by : 



p(z\R CRa i) oc ^-(z) x 
dz 



(1) 



dN, 



gal 



dM 



(M = i? cga i - (jl l (z) - K sal {z)) 



In this equation, dV c /dz is the derivative of the comov- 
ing volume with respect to z, and is the luminosity 
distance modulus. The cosmological parameters (Om; Oa) 
are set to the value (0.3; 0.). Taking the values (Om; Oa) = 
(0.3; 0.7) based on distant type Ia supernovae (Riess et al. 
1998) ; Perlmutter et al. 1999) decreases the computed rate 



by less than 5%. 

The distribution of the number of galaxies per unit 
absolute magnitude dN ga \/dM is the Schechter law: 



dM 



oc l0- - 4 x(« +1 ) x ( M - M *)exp (_io-0-4x(M-M,)) (2 ) 



measured by the LCRS (Lin et al. 1996) where QfLCRS 
-0.7 and M* floLC RS 
iVlOOkms^Mpc" 1 ). 



= -20.29 + 51og/i (h = 



The K-correction is simply given by K sa ,\(z) = 
2.5 log(l + z) as proposed by Lin et al. ( 199£| ). 

For a galaxy magnitude of i? C gai — 17, the mean and 
rms value of the ^-distribution thus obtained are respec- 
tively - 0.1 and - 0.05. 

The model was checked by comparing, for each galaxy 
with an LCRS redshift, the mean of p(z\R c ga j) with the 
measured LCRS redshift. Satisfactory agreement was ob- 
tained after adding 0.25 to the EROS total magnitudes 
(as measured with SExtractor) so as to align them with 
the isophotal LCRS magnitudes. 

5. The detection efficiency and the supernova rate 

The supernova rate 1Z in the rest frame is expressed in 
SNu, i.e. in supernovae per unit time per unit blue lumi- 
nosity ISNu = 1SN /10 10 L Qs /100yrs. The number of 
supernovae M expected to be detected is given by a sum 
over the observed galaxies i weighted by their blue lumi- 
nosities Lf. 



M = n x s, 

S = y^ Lj I ei(t, Zi) dt 



(3) 
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where 6j(t, Zi) is the efficiency to detect in galaxy % a type 
la supernova whose maximum occurs at time t in the su- 
pernova rest frame. The efficiency clearly depends on the 
galaxy redshift zj but also on other factors such as the 
galactic luminosity and surface-brightness profile. The in- 
tegral / Cidt is sometimes called the "control time" for the 
galaxy i. 

We now need to determine the sum S so the rate can 
be evaluated through Eq. ||. We do not know the lumi- 
nosities of most of the observed galaxies because most of 
their redshifts are not known. As discussed in the previous 
section, we do, however, know the probability p(z\R CS al) 
that the galaxy of magnitude i? cga i has a redshift z. We 
therefore replace the above expression for S with 



Detection efficiency (1) 



E 



dz p{z\R ci ) Lt(R ci ,z) 



i(t,z)dt (4) 



where the blue luminosity is computed using the mag- 
nitude R c i and the redshift z. The double integral in this 
equation was evaluated for each galaxy by Monte Carlo 
simulation of supernovae that are subjected to the same 
detection procedures as the real data. 

Because there is a better match between the EROS 
visible band and the standard V filter than with the stan- 
dard B filter, the absolute luminosities of galaxies are first 
computed in the standard V band, and then converted 
into absolute luminosities in the standard B band. We 
use the difference between the solar < B — V > colour: 
<B— V>q= 0.65 and a mean <B—V> colour for galaxies, 
compu ted from the results of Fioc & Rocca-Volmerange 
( |1997| ): <cB-^> gal = 0.77. 

5.1. Supernova simulation 

The Monte Carlo simulation proceeds as follows. For each 
potential host galaxy, a series of redshifts is drawn ran- 
domly according to the distribution appropriate for the 
galaxy apparent magnitude. For each redshift, a time t 
since maximum is drawn uniformly between —15 and +120 
days in the supernova rest frame. The absolute magnitudes 
My j and Mr c of the supernova at that time are taken from 
the templates in these standard filters provided by Riess 
et al. ( 1996j ). To account for the observed spread of su- 
pernova luminosities, a gaussian scatter AM = 0.2 mag 
was added. The apparent magnitude of the simulated 
supernova is computed from the redshift using the K- 
correction from Nugent & Kim (in preparation) and as- 
suming (S1m,Qa) = (0.3,0.) as above. 

The supernova standard apparent magnitudes are then 
converted into a Vetos flux using our calibration. The posi- 
tion of the supernova in the host galaxy is drawn according 
to a two-dimensional Gaussian distribution, the first and 
second order moments and orientation of which are those 
of the host object. 

As supernovas are stellar sources, they appear on the 
frame with the same PSF as the field stars do. The PSF 




22 23 
V magnitude 

Fig. 3: The detection efficiency e as a function of the simu- 
lated supernova magnitude Vj at discovery. It corresponds 
to a limiting magnitude Vj ~ 21.5. The maximum effi- 
ciency is below 1 because of the masked areas near CCD 
defects and saturated stars. 



Detection efficiency (2) 



6 0.35 F" 

0.3 '- 

0.25 

0.2 i- 

0.15 E- 

0.1 E- 

0.05 r 

E-l 



1 2 3 4 5 6 

dist. from host core (arcsec) 

Fig. 4: Dependence of the detection efficiency e on the dis- 
tance between the simulated supernova and its host galaxy 
core. The average efficiency is as low as ^30% because 
supernova light curves are sampled up to 120 days after 
maximum, causing ~ 75% of them to drop below Vj = 22. 



is modeled by a two-dimension integrated Gaussian. Its 
second order moments are set to the mean moments com- 
puted by fitting the stars on the same CCD quadrant, in 
order to take into account the spatial variation of the PSF. 

5.2. Detection efficiency 

The detection efficiency is derived by comparing the list 
of simulated supernovae and the list of the supernovae se- 
lected by the search cuts. 

The detection efficiency as a function of magnitude is 
shown in Fig. |3|. It appears as a smoothed step function, 
corresponding to a limiting magnitude Vj ~ 21.5. There 
is a slight dependence on the distance of the supernova 
from the host galaxy core, as shown in Fig. H. 
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Detection efficiency (3) 



E 
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0.8 
0.6 
OA 
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host surface brightness 
R./arcsec 2 




5 10 15 20 25 
host local relative gradient 



Fig. 5: Detection efficiency e as a function of the host 
galaxy surface brightness (left pannel) and of the local 
gradient at the point where the simulated supernova is 
added (right pannel). The local gradient was taken to be 
g = max(pixel(a; i ,y i ) -pixel(x J -,y J ))| j ._ i |< 1 /S where B 
is the local Poisson noise. The distributions of the sur- 
face brightness and the local gradient are superimposed 
(dashed lines). The distributions are shown for supernovae 
whose magnitude is in the range 20.5 < Vj < 21. 



As foreseen, the detection efficiency depends on the 
host galaxy characteristics, such as the surface brightness 
or the local gradient where the simulated supernova is 
added. These results are presented in Fig. bL where the 
distribution of these two parameters have been superim- 
posed on the efficiency. The distributions are shown for su- 
pernovae whose magnitude is in the range 20.5 < Vj < 21. 
They peak where the efficiency is good, thus moderating 
the impact of the efficiency behaviour. 



5.3. Redshift distribution 

The Monte-Carlo computed redshift distribution of de- 
tected supernovae is shown in Fig. |[ The thin solid line 
shows the distribution when all the analysis cuts are ap- 
plied. For comparison, the dotted line shows the redshift 
distribution when only two cuts are applied, one on the su- 
pernova magnitude Vj < Vj = 21.5, and one on the host 
galaxy V E ros magnitude equivalent to R c g a i 

gain = 

18.7. The thick solid line shows the distribution obtained 
with an analytic model where we assume that the type la 
supernova rate is proportionnal to the host galaxy lumi- 
nosity and that the galaxy luminosity distribution is given 
by the LCRS Schechter law (Eq. 2). In this analytic model 
the cuts are described by two step- functions whose thresh- 
olds are respectively Vj and i? C gai - Using this analytic 
model, we find that imposing the host galaxy magnitude 
cut reduces the number of detected supernova by a factor 
- 2.3. 



Expected redshift distribution 
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Fig. 6: The expected redshift distribution. The thin solid 
line shows the distribution in z derived from the Monte- 
Carlo simulation when the analysis cuts are applied. The 
dotted line shows the same but with only two cuts ap- 
plied, one on the supernova magnitude (Vj < 21.5) and 
the other one on the host galaxy Vei-os magnitude equiv- 
alent to R c ga i < 18.5. The latter distribution compares 
well with a simple analytic model as described in the text 
(thick solid line). The arrows indicate the redshift values 
of the 4 discovered type la supernovae. 



The observed redshift distribution of the 4 discov- 
ered type la supernovae is in agreement with the Monte- 
Carlo distribution (the Kolmogorov-Smirnov compatibil- 
ity test gives a 25% probability), although the small- 
number statistics prevents any definitive conclusion. The 
mean and rms of the expected redshifts, respectively 0.14 
and 0.06, are to be compared with the mean and rms of 
the observed redshifts, respectively 0.10 and 0.06. 



5.4- Explosion rate 

The rate is computed by dividing the number of observed 
supernovae J\f b s by the sum S : 



K = M ohs /S 



(5) 



The sum obtained is S = 9.09 10 12 /i~ 2 L 0B yr. For 
A/"obs = 4 type la supernovae, this corresponds to a rate 



K = 0.44+^° h 2 /10 1U L Q£ >/ 100 yr where the error is sta- 
tistical and given at a 68% confidence level. 



6. Discussion 

The numerical calculation of the efficiency makes it pos- 
sible to study various systematic effects. We study here 
the effects inherent in the hypotheses used in the Monte 
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Carlo simulation, such as the calibration relations or the 
assumed type la supernova light-curve shape. 

6.1. Calibration 

The calibration enters the calculation at three different 
steps: when measuring the galaxy magnitude R c , which 
is used to draw its redshift, when computing the ADU 
flux of a simulated supernova as a function of its Vj and 
R c magnitudes, and when estimating the galaxy absolute 
luminosity in the Vj filter for the rate computation in SNu. 
The redshift galaxy distribution is verified with the LCRS 
galaxies, therefore we shall only take into account here the 
last two steps. Their contributions are of opposite sign, 
so that the estimated 10% uncertainty in the calibration 
zero-point yields only a 5% uncertainty in the computed 
rate. 

6.2. Galaxy counts and redshift distribution 

To take into account the requirement for the detected su- 
pernova to be found in a galaxy, potential host galaxies 
were selected by requiring that their magnitude R c sat- 
isfy R c < R c0 (section 4.1). A 0.2 magnitude shift on R c0 
translates into a 10% shift on the rate. 

The galaxy sample misidentification is estimated to 
be about 10% (section 4.1), which gives rise to a 10% 
uncertainty in the computed rate. 

6.3. Supernova luminosity distribution 

In the simulation, the supernova light curve was approx- 
imated by Vj and R c templates of Riess et al. ( |1996| ). 
The light-curve shape is described through two parame- 
ters: the peak magnitude Mvj max and the stretch factor 
s (Perlmutter et al. 1999 ). We assumed in the analysis 
that Mv 7 max = —18.4 + 51og/i with an added gaussian 
scatter of 0.2 mag, which corresponds to the observed dis- 
tribution for the 16 supernovae from the Hamuy et al. 
( |l996| ) sample reanalysed in Perlmutter et al. ( 1999 ). No 
stretch factor correction was taken into account, i.e. s was 
set to s = 1. To test the robustness of this analysis, we 
reperformcd the analysis assuming, following Perlmutter 
et al. ( p-999 ), that the peak magnitude is described by 



Tabic 3: Summary of systematic errors. 



M Vj max = M Vj ma x - 0.6 X (s - 1 ) , where the M Vj max dis- 
tribution is a Gaussian centered on Myj = —18.4 + 5 log/i 
with a scatter a — 0.2 mag., and s is a Gaussian centered 
on s — 1 with a scatter a = 0.1. The rate calculated with 
these assumptions differed from that calculated with the 
standard assumptions by less than 1%. 

On the other hand, it seems possible that the type la 
supernova mean peak magnitude might be uncertain to 
about 0.1 mag. A shift of 0.1 mag in the peak magnitude 
results in a 5% variation of the rate. 



Magnitude calibration 


± 5% 


Limiting magnitude R c0 for the 
galaxy selection 


±10% 


Galaxy misidentification 


±10% 


SN mean peak magnitude 


± 5% 


SN type misidentification 


+25 o/ 
- % 


Total systematic error 


+30 <y 
-16 % 



6.4- Total number of detected type la 

The type of SN 1997eb is unclear. This is taken into ac- 
count by adding i ' o SNu to the systematic uncertainty. 



7. Conclusion 

We measure a rate of: 



n 44 +0.35 +0.13 1 2 
u -^-0.21 -0.07 11 



SNu 



where the systematic errors from Tab. || have been added 
in quadrature. The statistical error corresponds to a 68% 
confidence level. This rate compares well with the rate 
derived by Cappellaro et al. ( |1999| ) from the combined 
sample of five nearby photographic searches, TZ — 0.35 ± 
0.11 h 2 SNu at a mean redshift of z ~ 0.01. Similarly, our 
rate compares well with the Supernova Cosmology Project 
rate for type la at z ~ 0.4 derived by Pain et al. (1996), 
TZ = 0.82+°;^ h 2 SNu. 

The rate can also be expressed in /i 3 /Mpc 3 /yr by mul- 
tiplying the rate in SNu by the luminous density of the 
universe p L = f .4 ± 0.1 10 8 h L Mpc" 3 
This gives: 



(Lin ct al. |1996[ ) 



TZ = 0.62 ± 



+0.49+0.19 
-0.29-0.11 



10 -4 /i 3 /Mpc7yr. 
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